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Abstract. The aim of these notes is to explore possible ways of extending 
Smirnov's proof of Cardy's formula for critical site-percolation on the trian- 
gular lattice to other cases (such as bond-percolation on the square lattice); 
the main question we address is that of the choice of the lattice embedding 
into the plane which gives rise to conformal invariance in the scaling limit. 
Even though we were not able to produce a complete proof, we believe that 
the ideas presented here go in the right direction. 
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Introduction 

It is a strongly supported conjectured that many discrete models of random media, 
such as e.g. percolation and the Ising model, when taken in dimension 2 at their 
critical point, exhibit conformal invariance in the scaling limit. Indeed, the uni- 
versality principle implies that the asymptotic behavior of a critical system after 
rescaling should not depend on the specific details of the underlying lattice, and in 
particular it should be invariant under rotations (at least under suitable symmetry 
conditions on the underlying lattice). Since by construction a scaling limit is also 
invariant under rescaling, it is natural to expect conformal invariance, as the local 
behavior of a conformal map is the composition of a rotation and a rescaling. 

On the other hand, conformally invariant continuous models have been thor- 
oughly studied by physicists, using tools such as conformal field theories. In 2000, 
Oded Schramm ([H]) introduced a one-parameter family of continuous bidimen- 
sional random processes which he called SLE processes, as the only possible scaling 
limits in this situation, under the assumption of conformal invariance; connections 
between SLE and CFT are now quite well understood (see e.g. Oil]). 



However, actual convergence of discrete models to SLE in the scaling limit 
is known for only a few models. The case on which we focus in this paper is that 
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of percolation. The topic of conformal invariance for percolation has a long history 
— see [13] and references therein for an in-depth discussion of it. 

In the case of site-percolation on the triangular lattice, it is a celebrated result 
of Smirnov ([H>]) that indeed the limit exists and is conformally invariant. While 
the proof is quite simple and extremely elegant (see section [3] below and references 
therein), it is very specific to that particular lattice, to the point of being almost 
magical; it is a very natural question to ask how it can be generalized to other 
cases, and in particular to bond-percolation on the square lattice. Universality and 
conformal invariance have indeed been tested numerically for percolation in various 
geometries (see e.g. [12]), and conformal invariance (assuming the existence of the 
limit) is known in the case of Voronoi percolation (see [4]). 

In fact, it seems that the question of convergence itself has hardly been ad- 
dressed by physicists, at least in the CFT community — a continuous, conformally 
invariant object is usually the starting point of their work rather than its outcome. 
Techniques such as the renormalization group do give reason to expect the exis- 
tence of a scaling limit and of critical exponents, but they seem to not give much 
insight into the emergence of rotational invariance. 

This is not surprising in itself, for a very trivial reason: Take any discrete 
model for which you know that there is a conformally invariant scaling limit, say 
a simple random walk on Z 2 , and deform the underlying lattice, in a linear way, 
so as to change the aspect ratio of its faces. Then the scaling limit still exists 
(it is the image of the previous one by the same transformation); but obviously 
it is not rotationally invariant. Since all the rescaling techniques apply exactly 
the same way before and after deformation, they cannot be sufficient to derive 
rotational invariance. A trace of this appears in the most general statement of 
the universality hypothesis (see e.g. [131 section 2.4]): To paraphrase it, given any 
two periodic planar graphs, the scaling limits of critical percolation on them are 
conjugated by some linear map g. 

The main question we address in these notes is the following: Given a discrete 
model on a doubly periodic planar graph, how to embed this graph into the plane 
so as to make the scaling limit isotropic? If the graph has additional symmetry (as 
for instance in the case of the square or triangular lattices), the embedding has 
to preserve this symmetry; so a restating of the same question in the terms of the 
universality hypothesis would be, absent any additional symmetry for one of the 
two graphs involved, can one determine the map gl 

The most surprising thing (to me at least) about the question, besides the 
fact that it appears to actually be orthogonal to the interests of physicists in that 
domain, is that its answer turns out to depend on the model considered. In other 
words, there is no absolute notion of a "conformal embedding" of a general graph. 
In the case of the simple random walk, the answer is quite easy to obtain, though 
it does not seem to have appeared in the literature in the form we present it here; 
in the case of percolation, I could find no reference whatsoever, the closest being 
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the discussion and numerical study of striated models in [13] where, instead of 
looking at a different graph, the parameter p in the model is chosen to depend on 
the site in Z 2 in a periodic fashion — which admittedly is a very related question. 

The paper is roughly divided into two parts. In the first one, comprised of 
the firs two sections, we introduce some notation and the general framework of 
the approach, and we treat the case of the simple random walk. This is enough to 
prove that the correct embedding is not the same for it as for percolation; we then 
argue that circle packings might give a way to answer the question in the latter 
case. In the second part, which is of a more speculative nature, we investigate 
Smirnov's proof in some detail, and rephrase it in such a way that its general 
strategy can be applied to general triangulations. We then describe the two main 
steps of a strategy that could lead to its generalization, though we were able to 
perform none of the two. 

1. Notation and setup 
1.1. The graph 

We first define the class of triangulations of the plane we are interested in. Let T 
be a 3-regular finite graph of genus 1 (i.e., a graph that is embeddable in the torus 
T 2 := R 2 /Z 2 but not in the plane, and having only vertices of degree 3). For ease 
of notation, we assume that T is equipped with a fixed embedding in T 2 , which 
we also denote by T. The dual T* of T (which we also assume to be embedded in 
the torus once and for all) is then a triangulation of T 2 . 

Let T (resp. T*) be the universal cover of T (resp. T*): Then T and T* are 
mutually dual, infinite, locally finite planar graphs, on which Z 2 acts by translation. 
We are interested in natural ways of embedding T into the complex plane C. Let 
Ti (the meaning of the notation will become clear in a minute) be the embedding 
obtained by pulling T back using the canonical projection from R 2 to T 2 — we 
will call Ti the square embedding of T. 

For every a € C \ R, let ip a : C — > C be the R-linear map defined by 
(p a (x + iy) = x + ay (i.e., it sends 1 to itself and i to a) and let T a be the image 
of Ti by ip a . For lack of a better term, we will call T a the embedding of modulus a 
of T in the complex plane. 

Notice that the notation T a depends on the a priori choice of the embedding 
of T in the flat torus; but, up to rotation and scaling, the set of proper embeddings 
of T obtained starting from two different embeddings of T is the same, so no 
generality is lost (as far as our purpose in these notes is concerned). 

One very useful restriction on embeddings is the following: 

Definition 1. We say that an embedding T a of T in the complex plane is balanced 
if each of its vertices is the barycenter (with equal weights) of its neighbors; or, 
equivalently, if the simple random walk on it is a martingale. 
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Proposition 2. Let T be a i-regular graph of genus 1; Then, for every a £ H, there 
is a balanced embedding ofT in the complex plane with modulus a. Moreover, this 
embedding is unique up to translations of the plane. 

Proof. We only give a sketch of the proof, because expanding it to a full proof is 
both straightforward and tedious. The main remark is that any periodic embedding 
which minimizes the sum 52, over a period, of the squared lengths of its edges is 
balanced: Indeed, the gradient, with respect to the position of a given vertex, of 
S2 is exactly the difference between this point and the barycenter of its neighbors. 
(This would be true in any Euclidean space.) It is easy to use a compactness 
argument to prove the existence of such a minimizer. 

To prove uniqueness up to translation is a little trickier, but since it is not 
necessary for the rest of this paper, we allow ourselves to give an even sketchier 
argument. First, one can get rid of translations by assuming that a fixed vertex 
of T is put at the origin by the embedding; the set of all possible embeddings of 
modulus a is then parameterized by 2(|V r (T)| — 1) real-valued parameters, which 
are the coordinates of the locations of the other vertices in one period of T. In 
terms of these variables, S2 is polynomial of degree 2. It is bounded below by the 
squared length of the longest edge in the embedding, which itself is bounded below, 
up to a constant depending only on the combinatorics of the graph, by the square 
of the largest of the 2(|V(T)| — 1) parameters; so it goes to infinity uniformly at 
infinity. This implies that its Hessian (which is constant) is positive definite, so S2 
is strictly convex as a function of those variables. This immediately implies the 
uniqueness of the minimizer. □ 

An essential point is that, even though our proof uses Euclidean geometry, 
the fact that the embedding is balanced is a linear condition. In particular, if the 
embedding 1} is balanced, then so are all the other T a . The corresponding a priori 
embedding of T itself into the flat torus T 2 (which is also unique up to translations) 
will be freely referred to as the balanced embedding of T into the torus. 

1.2. The probabilistic model 

We will be interested in critical site-percolation on the triangulation T*; more 
specifically, the question we are interested is the following. Let S7 be a simply 
connected, smooth domain in the complex plane, and let A, B, C and D be four 
points on its boundary, in that order. For every 5 > 0, let fig be the largest 
connected component (in terms of graph connectivity) of the intersection of SI 
with ST a , and SI $ be its dual graph. fig should be seen as a discretization of SI at 
scale S. Let Ag, Bg, Cg and Dg be the vertices of fig that are closest to A, B, C 
and D respectively. 

The model we are most interested in is critical site-percolation on S1J; how- 
ever, most of the following considerations remain valid for other lattice models. 
Let Cg(Q, A, B, C, D) be the event that there is an open crossing in QJ, between 
the intervals AgBg and CgDg of its boundary. Under some symmetry conditions 
on T, Russo-Seymour- Welsh theory ensures that at criticality, the probability of 
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Cs(Cl, A, B, C, D) is bounded away from both and 1 as 8 goes to 0. Its limit was 
conjectured by Cardy (see [6]) using non-rigorous arguments from conformal field 
theory; actual convergence was proved, in the case of the triangular lattice (embed- 
ded in such a way that its faces are equilateral triangles), by Smirnov (see |15[l3]). 
We defer the statement of the convergence to a later time. The following definition 
has become standard: 

Definition 3. Assume that, for every choice of (0, A, B, C, D), the probability of 
the event Cs(Q, A, B, C, D) has a limit / a (^, A, B,C,D) as S — ► — we will refer 
to this by saying that the model has a scaling limit. We say that the model is 
conformally invariant in the scaling limit if, for every conformal map $ from fl to 
$(0), one has 

A, B, C, D) = / a (*(fi), $(A), <&(C), *(!>)). 

This is equivalent to saying that f a (Q,A,B,C,D) only depends on the modulus 
of the conformal rectangle (ft, A, B, C, D). 

(Notice that the extension of <I> to the boundary of Q, which is necessary 
for the above definition to make sense, is ensured as soon as Q, is assumed to be 
regular enough.) 

2. Periodic embeddings 
2.1. Uniqueness of the modulus 

Given T, it is natural to ask whether it is possible to choose a value for a which 
provides conformal invariance in the scaling limit. There are two possible strategies: 
Either give an explicit value for which "a miracle occurs" (in physical terms, for 
which the model is integrable — this is what Smirnov did in the case of the 
triangular lattice), or obtain its existence in a non-constructive way — which is 
what we are trying to do here. 

A reassuring fact is that, whenever such an a exists, it is essentially unique: 

Proposition 4. For every graph T , there are either zero or two values of a such 
that critical site-percolation on T* is conformally invariant in the scaling limit. In 
the latter case, the two values are complex conjugates of each other. 

Proof. The key remark is the following: Let (3 be a non-real complex number. Since 
the event Cs is defined using purely combinatorial features, one can push the whole 
picture forward through (pp without changing its probability. Let a' = (pp(a): (pp 
then transforms f2 into ipp{£l) and the lattice T a into T a >. So, assuming convergence 
on both sides, one always has 

/ a (n, A, B, C, D) = fa>(<Pp(Q>), ¥>/3{A), MB), V p{C),ip {D)). 

In the case /3 = —i, <pp is simply the map z i— ► z. In that case, the modulus of 
the conformal rectangle (ip-i{£l),D, (7, -B, A) is the same as that of (O, A, B, C, D), 
and clearly the event Cs is invariant when the order of the corners is reversed. So, 
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conformal invariance for T a and the previous remark implies that /a(fl, A, B, C, D) 
still only depends on the modulus of the conformal rectangle — in other words, 
if critical percolation T a is conformally invariant in the scaling limit, that is also 
the case on T a . 

Now assume conformal invariance in the scaling limit for two choices of the 
modulus in the upper-half plane; these moduli can always be written as a and 
a' = tpp(a) for an appropriate choice of (3 € H\ {i}. Still using the above remark, 
all that is needed to arrive to a contradiction is to show that f a does actually 
depend on the modulus of the rectangle (i.e., that it is not constant), and that 
there exist two conformal rectangles with the same modulus and whose images by 
ipp have different moduli. 

For the former point, it is enough to prove that for every choice of p, p' > 0, 
the probability of crossing the rectangle [0, p] x [0, 1] horizontally is strictly larger 
than that of crossing [0, p + p'] x [0, 1]. This is obvious by Russo-Seymour- Welsh: 
The event that there is a vertical dual crossing in ST* n [p + S, p + p'\ x [0, 1] is 
independent of Cs{[0, p] x [0, 1], p, p + i, i, 0) and its probability is bounded below, 
uniformly in S < p'/lO, by some positive e depending only on p and p' . Hence, 
still assuming that the limits all exist as 5 — > 0, 

f a ([0,p + p'\ x [0,l],p + p',p + p' +i,i,0) (l-e)/ a ([0,p] x [0,l],p,p + i,i,0). 

For the latter point, assume that •pp preserves the equality of moduli of con- 
formal rectangles. Let Q = [0, l] 2 be the unit square. By symmetry, the conformal 
rectangles (Q, 0, 1, 1 + i, i) and (Q, 1, 1 + i, i, 0) have the same modulus; on the 
other hand <pp(Q) is a parallelogram, and by our hypothesis on (pp it has the same 
modulus in both directions. This easily implies that it is in fact a rhombus. If now 
Q' is the square with vertices 1/2, 1 +i/2, 1/2 + i, i/2, <pp{Q') is both a rhombus 
(by the same argument) and a rectangle (because its vertices are the midpoints of 
the edges of <pp(Q) which is a rhombus). Hence pp(Q') is a square, and so is <pp{Q), 
and in particular (3 = <pp(i) = i, which is in contradiction with our hypothesis. □ 

When such a pair of moduli exists, we will denote by a^ erc the one with 
positive imaginary part. The same reasoning can be done for various models, and 
in each case where the scaling limit exists and is non-trivial, there will be a pair of 
moduli making it conformally invariant; we will distinguish them from each other 
by using the name of the model as a superscript (so that for instance aj. w makes 
the simple random walk conformally invariant in the scaling limit — cf. below). 

When an argument does not depend on the specific model (as is the case in 
the next subsection), we will use the generic notation or as a placeholder. 

2.2. Obtaining ar by symmetry arguments 

It should be noted that, because the value of ar (when it exists) is uniquely 
defined by the combinatorics of T, there are cases where additional symmetry 
specifies its value uniquely. Indeed, assume ^ is a graph isomorphism of T which 
is neither a translation nor a central symmetry; for every a, it induces a topological 
isomorphism of T a . Assume without loss of generality that the origin of the plane 
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is chosen to be one of the vertices of T a ; let z = ^(0), z\ — \&(1) and z a = ^(a) 
(notice that both 1 and a are also vertices of T a ). 

Assume a = a^ crc . Because 4" is an isomorphism, it preserves site-percolation; 
so, in particular, critical site-percolation on ^{T a ) is conformally invariant in the 
scaling limit. By Proposition [U this implies that 

*(a) - *(0) 



G {a, a}. 



(2.1) 



21 - z (1) - *(0) 

This condition is then enough to obtain the value of or- There are two natural 
examples of that (illustrated in Figures [TJ and [2]) , which we now describe. 





Figure 1. The graphs T h (left) and T s (right), embedded into 
T 2 in a balanced way with a vertex at the origin; empty circles 
and dotted lines represent the dual graphs. Both are represented 
using their square embedding, so the triangles in T/, are not equi- 
lateral. 




Figure 2. The same graphs as in Figure [J, with the origin on a 
vertex of the dual. 



• Let Th be one period of the honeycomb lattice, embedded into T 2 in such 
a way that every vertex is the barycenter of its neighbors (we will call such 
an embedding balanced); since we take T 2 to be a square, the coordinates of 
the vertices of T h are (0,0), (1/3,0), (1/2,1/2) and (5/6,1/2). There is an 
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isomorphism of order 3 of f^, corresponding to rotation around (0, 0); on 
T a , it sends to z = 0, 1 to Z\ = (3a — l)/2 and a to z a = —(1 + a)/2. 
Since preserves orientation, Equation l|2.ip leads to 



the faces of T f * as equilateral triangles, and those of Th as regular hexagons. 
• Let T s be chosen in such a way that T* has the topology of the centered square 
lattice; if again the embedding is balanced, the coordinates of the vertices of 
T s are (0,0), (1/2,0), (1/4,1/4) and (1/4,3/4). There is an isomorphism 
4" s of order 4 of T s , corresponding to a rotation rotation around the vertex 
(1/4,0) of T*. In that case 



so ar s — i- Again not surprisingly, this corresponds to the usual embedding 
of the square lattice using — well, squares. 

Of course, identifying ar in those cases is a long way from a proof of conformal 
invariance; but it would seem that understanding, in the general case, what a^ Grc 
is would be a significant progress in our understanding of the process. 

2.3. Embedding using random walks 

As an aside, in this subsection and the next we describe two natural ways of embed- 
ding a doubly periodic graph into the complex plane, which both have something 
to do with conformal invariance. 

Let T be a finite 3-regular graph of genus 1, embedded in T 2 in a balanced 
way, and let (X n ) n ^o be a simple random walk on it. For simplicity, assume that 
(X n ) is irreducible as a Markov chain. (Both 3-regularity and irreducibility are 
completely unnecessary as far as the results presented here are concerned, and 
the same reasoning would work in the general case, but notation would be a little 
tedious.) Since T is finite, (X n ) converges in distribution to the unique invariant 
measure, which, because T is 3-regular, is the uniform measure on V(T); moreover 
the convergence is exponentially fast. 

Now pick a£l, and lift (X n ) to a simple random walk (Z n ) on T a . By the 
balance condition on the embedding, it is easy to check that (Z n ) is a martingale; 
exponential decay of correlations between its increments is enough to obtain a 
central limit theorem (cf. for instance [9] and references therein). To write the 
covariance matrix in a convenient form, we need some notation. For each (oriented) 
edge e of T, choose zi(e) and 22(e) in T a in such a way that they are neighbors 
and the edge (z 1 (e) 7 z 2 (e)) is a pre-image of e by the natural projection from T a 




z a - z = (-3/4 - q/4) - (1/4 - a/4) = -1 
zi-zo (l/4 + 3a/4) - (1/4 -a/4) a 
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to T; let e a := z 2 (e) — Zi(e) — obviously it does not depend on the choice of Z\(e) 
and 22(e). Define 

1 



£™(T) := 
E**(T) := 



\E(T)\ 
1 

1 



E (3?e Q ) 2 , 

e£B(T) 

E ( 3e -) 2 ' 

ee£(T) 

E (Ke«)(9e a ). 



eGE(T) 

It is not difficult to compute the covariance matrix of the scaling limit of the walk: 

Proposition 5. As n goes to infinity, n~ x l 2 Z n converges in distribution to a Gauss- 
ian variable with covariance matrix 

- Ys x a x {T) E2*(T) " 
ES"(T) S^(T) 



S«(T) 



Proof. The walk is centered by definition; the existence of a Gaussian limit is a 
direct consequence of the exponential decay of step correlations. All that remains 
to be done is to compute the covariance matrix. We focus on the first matrix entry, 
the others being similar. We have 



Var(n _1 / a RZ n ) = E 



1 n—1 

-j2E\mz k+1 -z k )y 



(the other terms disappear by the martingale property). We know that Z k +i — Z k 
converges in distribution, because the walk on T converges in distribution; its limit 
is the distribution of e a where e is an edge of T chosen uniformly. By Cesaro's 
Lemma, the expression above then converges to S™; the computation of the other 
entries in H a (T) it exactly similar. □ 

Even though the previous definition of conformal invariance in the scaling 
limit does not apply directly in this case, its natural counterpart is to ask for 
the scaling limit of the walk to be rotationally invariant (i.e., to be standard 
two-dimensional Brownian motion); this is equivalent to saying that the covari- 
ance matrix £ a (T) is scalar, and since its entries are real, yet another equivalent 
formulation is 



[ Sr(r) _ E w( T )] + ^ ( T)] = 



E 

e£E(T) 



0. 



The last equation is a second-degree equation in a with real-valued coeffi- 
cients. If q £ t, all the terms are non-negative and at least one is positive, so 
the equation has no solution in K; letting a go to +00 along the real line leads 
to |-B(T)| positive terms, at least one of which is of order a 2 , so the coefficient 
in a 2 in the equation is not zero. Hence the equation has exactly two solutions 
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which are complex conjugate of each other — the situation is very similar to the 
one in Proposition [H For further reference, we let a^ w be the one with positive 
imaginary part. One advantage of this choice (besides the fact that it exists for 
every doubly periodic graph) is that the value of is very easy to compute. 

Remark 6. In the more general case of a doubly periodic graph but without the 
assumptions of 3-regularity and irreducibility (but still assuming that the embed- 
ding is balanced) , the condition e a ls still necessary and sufficient for the walk 
to be isotropic in the scaling limit — and the proof is essentially the same, so we 
do not delve into more detail. 

Remark 7. Of course, in the cases where T has some additional symmetry, a^ w 
is the same as that obtained in the previous subsection using symmetry . . . 

Remark 8. One can also look at a simple random walk on the dual graph T*, and 
ask for which values of a this dual walk is isotropic in the scaling limit. As it turns 
out, the modulus one obtains this way is the same as on the initial graph, in other 
words 

RW _ RW 

This is a very weak version of universality, and unfortunately there doesn't seem 
to be a purely discrete proof of it — say, using a coupling of the two walks. 



There is another natural way to obtain the same condition. We are planning 
on studying convergence of discrete objects to conformally invariant limits, so it 
is a good idea to look for discrete-harmonic functions on T a (with respect to the 
natural Laplacian, which is the same as the generator of the simple random walk 
on T a ). The condition of balanced embedding is exactly equivalent to saying that 
the identity map is harmonic on T a ; it is a linear condition, so it does not depend 
on the value of a. 

The main difficulty when looking at discrete holomorphic maps is that the 
product of two such maps is not holomorphic in general. But we are interested in 
scaling limits, so maybe imposing that such a product is in fact "almost discrete 
holomorphic" (in the sense that it satisfies the Cauchy-Riemann equations up to 
an error term which vanishes in the scaling limit) would be sufficient. 

Whether the previous paragraph makes sense or not — let us investigate 
whether the map £ : z t-y z 2 is discrete-harmonic. For every z e T , we can write 

z'~z e£E z (T) e£E z {T) 

(the term in ^ ze a vanishes because the embedding is balanced). So, if Q is discrete- 
harmonic, summing the above relation over z £ T gives the very same condition 
X] e a = as before; in other words, a^ w is the embedding for which z t-y z 2 is 
discrete-harmonic on average. 
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As a last remark, let us investigate how strong the condition of exact har- 
monicity of C is; so assume that a is chosen in such a way that AC is identically 
0. Let e be any oriented edge of T; let e' := r.e and e" :— r 2 .e be the two other 
edges sharing the same source as e. The balance condition on the embedding plus 
harmonicity of £ imply the following system: 



Up to rotation and scaling, one can always assume that e a = 1, so the system 
reduces to e' a + e" a = —1 and (e' Q ) 2 + (e^) 2 = —1. Squaring the first of these two 
relations and substracting the second, one obtains e' a e'^ = 1, so e' a and e" a are the 
two solutions of the equation 

X 2 +X + l = () 
which implies that {e' Q ,e„} = {e ±27 ™/ 3 }. To sum it up: 

Proposition 9. The only 3-regular graph on which the map £ : z i— > z 2 is discrete- 
harmonic is the honeycomb lattice, embedded in such a way that its faces are regular 
hexagons. 

So, imposing £ to be harmonic not only determines the embedding, it also 
restricts T to essentially one graph; but in terms of scaling limits, the condition 
that £ is harmonic on the average makes as much sense as the exact condition. 

2.4. Embedding using circle packings 

There is another way to specify essentially unique embeddings of triangulations, 
which is very strongly related to conformal geometry, using the theory of circle 
packings. It is a fascinating subject in itself and a detailed treatment would be 
outside of the purpose of these notes, so the interested reader is advised to consult 
the book of Stephenson [18] and the references therein for the proofs of the claims 
in this subsection and much more. 

We first give a version of a theorem of Kobe, Andreev and Thurston, special- 
ized to our case. It is a statement about triangulations, which is why we actually 
apply it to T* instead of directly to T. Notice that we do not assume T to be 
already embedded into the torus T 2 . 

Theorem 10 (Discrete uniformization theorem [TBI P- 51]). Let T* be a finite 
triangulation of the torus, and let T* be its universal cover. There exists a locally 
finite family (C v ) v£ v{f*) °f disks of positive radii and disjoint interiors, satisfying 
the following compatibility condition: C v and C v > are tangent if, and only if, v and 
v' are neighbors in T* . 

Such a family is called a circle packing associated to the graph T* . It is es- 
sentially unique, in the following sense: If (C' v ) is another circle packing associated 
to T* , then there is a map <p : C — ► C, either of the form z i— > az + b or of the 
form z t-^ az + b, such that for every v G V(T*), C' v = (p(C v ). 









(2.2) 
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Remark 11. The "existence" part of the above theorem remains true in a much 
broader class of graphs; essentially all that is necessary is bounded degree and 
recurrence of the simple random walk on it. (One can see that a packing exists by 
completing the graph into a triangulation.) The "uniqueness" part however fails in 
general, as is made clear as soon as one tries to construct a circle packing associated 
to the square lattice . . . 

A consequence of the uniqueness part of the theorem is the following: Let 
9 : T* — > T* be a translation along one of the periods of T* , and let C' v := Cg(„); 
according to the theorem, let <p be such that C' v — tp(C v ) for all v. Up to composition 
of ip by itself, one can always assume that it is of the form (p(z) = az + b. By the 
assumption of local finiteness of the circle packing, one has \a\ = 1; besides, the 
orbits of 9 are unbounded, so those of tp are too, and in particular it does not 
have a fixed point, which implies that a = 1 and b ^ 0. In other words, ip is a 
translation, i.e. the circle packing associated to T* is itself doubly periodic. 

As soon as one is given a circle packing associated to a planar graph, it 
comes with a natural embedding: Every vertex v € V(T*) will be represented by 
the center of C v , and if v' is a neighbor of V, the edge (v, v') will be embedded as a 
segment — which is the union of a radius of C v and a radius of C v > , because those 
two disks are tangent. One can then specify an embedding of T by putting each of 
its vertex at the center of the disk inscribed in the corresponding triangular face 
of (the embedding of) T*; the collection of all those inscribed disks is in fact a 
circle packing associated with the graph T (see Figure [3]). 




Figure 3. The circle packings associated to the graph T s (solid 
lines) and its dual (dotted lines). 
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Of more interest to us is the fact that the embedding is itself doubly peri- 
odic, by the previous remarks. Up to rotation, and scaling and maybe complex 
conjugation, one can assume that the period corresponding to the translation by 
1 (resp. i) in T 2 is equal to 1 (resp. a € H). Once again, the value of the modulus 
a is uniquely determined; for further reference, we will denote it by . 

Yet again, as soon as one additional symmetry is present in T, the value 
of a T v is the same as that obtained using the symmetry; this is again a direct 
consequence of the essential uniqueness of the circle packing. 

2.5. "Exotic" embeddings 

Looking closely at Smirnov's proof, one notices that essentially the only place 
where the specifics of the graph are used is in the proof of "integrability" or exact 
cancellation; we will come back to this in the next section, let us just mention 
that the key ingredient in the phenomenon can be seen to be the fact that ip(e) 
(as introduced earlier) is identically 0. This is equivalent to saying that all the 
triangles of the triangular lattice are equilateral. 

A way to try and generalize the proof is to demand that all the faces of T* 
be equilateral triangles. Of course this cannot be done by embedding it in the 
plane, even locally — the total angle around a vertex would be equal to 2ir only 
if the degree of the vertex is 6. But one can build a 2-dimensional manifold Mt 
with conic singularities by gluing together equilateral triangles according to the 
combinatorics of T*; since the average degree of a vertex of T* is equal to 6, the 
average curvature of the manifold (defined e.g. as the limit of the normalized total 
curvature in large discs) is 0. 

The manifold Mt is not flat in general (the only case where it is being 
the triangular lattice), but it is homeomorphic to the complex plane, and one 
can hope to see it as a perturbation of it on which some of the standard tools 
of complex analysis could have counterparts — the optimal being to be able to 
perform Smirnov's proof within it. This is no easy task, and is probably not doable 
anyway. 

To relate Mt to the topic of this section, one can try to define a module out 
of it. A good candidate for that is the following: Assume that Mt can be realized 
as a sub-manifold of R 3 (or in R d for d > 2 large enough) , in such a way that the 
(combinatorial) translations on T act by global translations of the ambient space, 
thus forming a periodic sub-manifold. Then there is a copy of 1? acting on it, and 
the afHne plane containing a given point of Mt and spanned by the directions of 
the two generators of that group is at finite Hausdorff distance from it; in other 
words, this realization of Mt looks like a bounded perturbation of a Euclidean 
plane. 

One can then look at the orthogonal projections of the vertices of T* (seen as 
points of Mt) onto that plane; this creates a doubly periodic, locally finite family 
of points of the Euclidean plane. It is not always possible to form an embedding 
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of f* in the plane from it (with disjoint edges); but it does define a value of a as 
above. 

Unfortunately, there are cases when this value of a is not well-defined, in the 
sense that it depends on the choice of Mt) this happens if the (infinite) polyhedron 
associated to T* , with equilateral faces, is flexible. The simplest example of this 
phenomenon is to take T to be two periods of the honeycomb lattice in each 
direction. 

2.6. Comparing different methods of embedding 

We now have at least two (forgetting about the last one) ways of giving a conformal 
structure to a torus equipped with a triangulation — which is but another way of 
referring to the choice of a. Assuming that critical percolation does have a scaling 
limit, it leads to a third choice o^ orc of it. 

It would be a natural intuition that all these moduli are the same, and corre- 
spond to a notion of conformal embedding of a triangulation (or a 3-regular graph) 
in the plane; and they all have a claim to that name. But this is not true in general: 
We detail the construction of a counterexample. Start with the graph T s and its 
dual T* ; and refine one of the "vertical" triangular faces of T* by adding a vertex 
in the interior of it, connected to its three vertices. In terms of the primal graph, 
this correspond to replacing one of its vertices by a triangle — see Figured! Let T' s 
be the graph obtained that way; we will refer to such a splitting as a refinement, 
and to the added vertex as a new vertex. 




Figure 4. Square (but not balanced) embeddings of T' s (solid) 
and its dual (dotted) ; the origin is taken as a point of T' s on the 
left, and as a point of the dual on the right, corresponding to the 
ones chosen for Figures [TJ and [2l 

In terms of circle packings, this changes essentially nothing; the new vertex of 
(T' s )* can be realized as a new disc without modifying the rest of the configuration 
(cf. Figured]). In terms of random walks, however, adding edges will modify the 
covariance matrix in the central limit theorem. The computation can be done 
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easily, as explained above, and one gets the following values: 




Figure 5. Splitting a face of a triangulation into 3 triangles, 
and the corresponding modification of its circle packing (added 
features represented by dashed lines) 



So, a RW and a CP are different in general. Is a Pcrc (provided it exists) one 
of them? An easy fact to notice is the following: Let T* be a triangulation of the 
torus and let (T")* be obtained from it by splitting a triangle into 3 as in the 
construction of T' s . Then, consider two realizations of site-percolation at p c = 1/2 
on both universal covers, coupled in such a way that the common vertices are in the 
same state for both models. In other words, start with a realization of percolation 
on T* and without changing site states, refine a periodic family of triangles of it 
into 3, choosing the state of each new vertex independently of the others and of 
the configuration on T* . 

If there is a chain of open vertices in T* , this chain is also a chain of open ver- 
tices in the refined graph — because all the edges are preserved in the refinement. 
Conversely, starting from a chain of open vertices in the refinement and removing 
each occurrence of a new vertex on it, one obtains a chain of open vertices in T*; 
the reason for that being that the triangle is a complete graph. Another way of 
stating the same fact is to say that opening (resp. closing) one of the new vertices 
cannot join two previously disjoint open clusters (resp. split a cluster into two 
disjoint components); they cannot be pivotal for a crossing event. 

Hence, the probability that a large conformal rectangle is crossed is the same 
in both cases (at least if the choice of the discrete approximation of its boundary 
is the same for both graphs, which in particular implies that it contains no new 
vertex), and so is f a for every choice of a (still assuming that it exists, of course). 
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If one is conformally invariant in the scaling limit, the other also has to be. In 
short. 



„ Perc rt Pcrc 
Otrp = a T , . 

Looking at circle packings instead of percolation, we get the same identity 
(as was mentioned in the particular case of T s ), with a very similar proof: Adding 
a vertex does not change anything to the rest of the picture, and we readily obtain 

CP CP 

a T = a T i . 

This leads us to the following hope, which we state as a conjecture even though it 
is much closer to being wishful thinking: 

Conjecture 12. Let T* be a triangulation of the torus. Then, the critical parameter 
{or site-percolation on its universal cover T* is equal to 1/2, and for every a£l, 
critical site-percolation on T* has a scaling limit. The value of the modulus a for 
which the model is conformally invariant in the scaling limit is that obtained from 
the circle packing associated to T* : 

„,Perc _ „,CP 



3. Critical percolation on the triangular lattice 

For reference, and as a way of introducing our general strategy, we give in this 
section a very shortened version of Smirnov's proof of the existence and conformal 
invariance of the scaling limit for critical site-percolation on the triangular lattice 
Th- The interested reader is advised to consult our previous note [3] for an "ex- 
tended shortening", or Smirnov's article [16] for the original proof; see the book of 
Bollobas and Riordan [5] for a more detailed treatment. Up to cosmetic changes, 
we follow the notation of [3]. 

Remark 13. Up to the last paragraph of the section, we are not assuming that 
the lattice we are working with is the honeycomb lattice; our only assumption is 
that we have an a priori bound for crossing probabilities of large rectangles which 
depends on their aspect ratio but not on their size (we "assume Russo-Seymour- 
Welsh conditions"). It is not actually clear how general those are; all the standard 
proofs require at least some symmetry in the lattice in addition to periodicity, but 
it is a natural conjecture that periodicity is enough. 

Here and in the remainder of this paper, r := e 2 ™/ 3 will be the third root 
of unity with positive imaginary part. Let T be a finite graph of genus 1, T a an 
embedding of modulus a of T in the complex plane; let V{T a ) (resp. E(T a )) be 
the set of vertices (resp. oriented edges) of T a . Each vertex z £ V(T a ) has three 
neighbours; let E z (T a ) be the set of the three oriented edges in E(T a ) having their 
source at z. That set can be cyclically ordered counterclockwise; if e € E z (T a ) is 
one of the three edges starting at z, we will denote by r.e (resp. r 2 .e) the next 
(resp. second to next) edge in the ordering. 
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Remark 14. In the particular case of the honeycomb lattice, seeing each edge as 
a complex number (being the difference between its target and its source), the 
notation r.e corresponds to complex multiplication by e 2 ^' 1 ^ 3 — in other words, 
r.e = re as a product of complex numbers. That is of course not the case in 
general, but we keep the formal notation for clarity. In what follows, whenever an 
algebraic expression involves the product of a complex number by an edge of T a or 
T*, as above the edge will be understood as the difference, as a complex number, 
between its target and its source; we will never use formal linear combinations of 
edges. The notation r.e (with a dot) will be reserved for the "topological" rotation 
within E z (T a ). 

Let again U be a smooth Jordan domain in the complex plane, and let A, 
B, C and D be three points on its boundary, in that order when following dVl 
counterclockwise. Let fls be the largest connected component of ft n ST a , and let 
As (resp. B$, C$, D$) be the point of Q,$ that is closest to A (resp. B, C, D). The 
main result in Smirnov's paper ([E]) is the following: 

Theorem 15 (Smirnov). In the case where T a is the honeycomb lattice, embedded 
so as to make its faces regular hexagons (i.e., when a = iy/3/3), critical site- 
percolation has a conformally invariant scaling limit. Iffl is an equilateral triangle 
with vertices A, B and C, then 



Knowing this particular family of values of f a is enough, together with con- 
formal invariance, to compute it for any conformal rectangle. The formula obtained 
for a rectangle is known as Cardy 's formula. 

To each edge e € E(T a ) corresponds its dual oriented edge e* G E(T*), 
oriented in such a way that the angle (e, e*) is in (0,vr). If — e denotes the edge 
with the same endpoints as e but the reverse orientation, then we have e** = — e. 
Define 



(where as above we interpret the edges e*, (r.e)* and (r 2 .e)* as complex numbers). 
It is easy to check that ip(e) = if, and only if, the face of T* corresponding to 
the source of e is an equilateral triangle; so, -0(e) can be seen as a measure of the 
local deviation between T a and the honeycomb lattice. An identity which will be 
useful later is the following: 



f a (Q,A,B,C,D) 



\CD\ 



\CA\- 



V>(e) := e* + r(r.e)* + r 2 (r 2 .e)* 





(3.1) 



e£E z (T a ) 



For every z & Slg, let Ea,s(z) be the event that there is a simple path of open 
vertices of Slg, joining two points of the boundary of the domain, which separates 
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z and A from B and C; let Ha '■= P[Ea,s(z)]- Define similar events for points B 
and C by a circular permutation of the letters, and let 



It is a direct consequence of Russo-Seymour- Welsh estimates that these functions 
are all Holder with some universal positive exponent, with a norm which does not 
depend on 5, so by Ascoli's theorem they form a relatively compact family, and as 
S — > they have subsequential limits which are Holder maps from Q to C; all that 
is needed is prove that only one such limit is possible. 

The key argument is to show that if h (resp. s) is any subsequential limit of 
(H$) (resp. (Ss)) as S — > 0, then h and s are holomorphic; indeed, assume for a 
moment that they are. Since s is also real-valued, it has to be constant, and its 
value is 1 by boundary conditions (e.g. at point A). On the other hand, along the 
boundary arc (A$B$) of dQg, Hc,s is identically 0, so the image of the arc (AB) 
by h is contained in the segment [l,r] of C; and similar statements hold mutatis 
mutandis for the arcs (BC) and (CA). By basic index theory, this implies that h 
is the unique conformal map sending Q to the (equilateral) triangle of vertices 1, r 
and t 2 , and that is enough to characterize it and to finish the proof of Theorem[T5l 

So, the crux of the matter, as expected, is to prove that the map h has to 
be holomorphic. The most convenient way to do that is to use Morera's theorem, 
which states that h is indeed holomorphic on ft if, and only if, its integral along 
any closed, smooth curve contained in ft is equal to 0. 

Let 7 be such a curve, and let 75 = (zo, Z\, ■ ■ ■ , zl s = zq) be a closed chain of 
vertices of fig which approximates it within Hausdorff distance 5 and has 0(5^ 1 ) 
points. Because the functions Hg are uniformly Holder, it follows that 



We want to prove that, for a suitable choice of a, the discrete integral on the 
left-hand side of that equation vanishes in the scaling limit. 

If e = (z, z') is an oriented edge of fig, define Pa,s{s) '■= P[EA,s(z')\EA,g(z)]; 
define Pb,s and Pc,g similarly. A very clever remark due to Smirnov, which is 
actually the only place in his proof where specifics of the model (as opposed to 
the lattice) are used, is that one can use color-swapping arguments to prove that, 
for every oriented edge, 



Sg(z) := H A ,g{z) + H B<5 (z) + H c ,g{z), 
Hg(z) := H A ,g(z) + tH b ^(z) + T 2 H c .g(z). 




PA,s(e) = P B ,s{T.e) = P c ,s(T 2 .e). 



(3-2) 



On the other hand, since differences of values of Hg between points of fig can 
be computed in terms of these functions P. : g, the discrete integral above can be 
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rewritten using them: Letting E(js) be the set of edges contained in the domain 
surrounded by 75 and using (|3.2p . one gets 



75 e&E{ ls ) 

A similar computation, together with the fact that e* + (r.e)* + (r 2 .e)* is identically 
equal to 0, leads to 



We again refer the reader to [3] for the details of this construction. 

Notice that it already implies that s is holomorphic, hence constant equal to 
1, independently of the value of a; so, whether h is holomorphic or not, it will send 
0. to the triangle of vertices 1, r and r 2 anyway. In the case of the triangular lattice 
embedded in the usual way, tfj(e) is also identically equal to 0, as was mentioned 
above, so h is itself holomorphic, and the proof is complete. 

The remainder of these notes is devoted to some ideas about how to extend 
the general framework of the proof to more general cases; it is not clear how close 
one is to a proof, but it is likely that at least one fundamentally new idea will 
be required. However, we do believe that the overall strategy which we will now 
describe is the right angle of attack of the problem. Do not expect to find any 
formal proof in what follows, though. 

4. Other triangulations 

4.1. Using local shifts 

The first natural idea when trying to generalize the construction of Smirnov is to 
try an apply it to more general periodic triangulations of the plane. Indeed, in all 
that precedes, up to and including Equation (|3.3p . nothing is specific to the regular 
triangular lattice, only Russo-Seymour- Welsh conditions (and their corollary that 
p c = 1/2) are needed. It is only at the very last step, noticing that tp was identically 
equal to 0, that the precise geometry was needed. 

The key fact that makes hope possible is the following (and it is actually 
similar to one of the points we made earlier): In the expression of the discrete 
integral as a sum over interior edges, each term is the product of two contributions: 

• ip( e ) which depends on the geometry of the embedding, and through that on 
the value of a; 

• Pa,s{b) which is only a function of the combinatorics of Qs- 

Even though fig as a graph does depend on the choice of a, one can make the 
following remark: Applying the transformation <pp (for some (3 £ H) to both the 
domain Q and the lattice ST a does not change Q a as a graph. In particular it does 
not change the value of Pa,s{&)- 




(3.3) 




(3.4) 
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One can then see the whole sum as a function j3, say Sn,s(P)- Because <pp{z) 
is a real-affine function of j3, so is Sn,s\ one can then try to solve the equation 
Sn,s(P) = in (3. Using the corresponding ipp, one gets a joint choice of a domain, 
a lattice modulus and mesh, and a curve 7 making the discrete contour integral 
vanish. 

It the modulus thus obtained actually did not depend on fl, S or 7, we would 
be done — call it cJp elc and there is only bookkeeping left to do. However we do not 
even know whether it has a limit as 6 I . . . An alternative is as follows. Because 
the lattice is periodic, it makes sense to first look at the sum ^{ e )PA,s{ e ) over 
one fundamental domain. If that is small, then over the copy of the fundamental 
domain immediately to the right of the previous one, the terms ip(e) are exactly 
the same, and one is lead to compare Pa,s for two neighboring pre-images of a 
given edge of T. 

So, let e be an edge of fls, and let e + 5 be its image by a translation of one 
period to the right. Making the dependency on the shape of the domain explicit in 
the notation, one can replace the translation of e by a translation of the domain 
itself and the boundary points in the opposite direction, to obtain 

P£ 5 (e + 5)=P£/(e). (4.1) 

To estimate the difference between this term and the corresponding one in fl, one 
can consider coupling two realizations of percolation, one on £ls and the other in 
tls — S, so that they coincide on the intersection between the two. 

The event corresponding to Pas( c ) ls that there is an open simple path 
separating the target of e and A from B, and C, and that no open simple path 
separates the source of e and A from B and C; this is equivalent to the existence 
of 3 disjoint paths from the 3 vertices of the face at the source of e to the 3 "sides" 
of the conformal triangle (f2, A, B, C), two of them being formed of open vertices 
and the third being formed of closed vertices — cf. Figure [6j For this to happen 
in fl but not in fl — S, one of these arms needs to go up to dfl but not to d(fl — S), 
and the only way for this to be realized is for a path of the opposite color to 
prevent it; this can be done in finitely many ways, Figure [6] being one of them; 
Pas ~ Pas S can then be written as the linear combination of the probabilities 
of finitely many terms of that form — half of these actually corresponding to the 
reversed situation, where arms go up to d(Q — 5) but not up to dtt. 

In the case corresponding to Figure El and all the similar ones, one sees that 
3 arms connect the source of e to the boundary of ft n (SI — 5), and on at least one 
point of that boundaries there have to be 3 disjoint arms of diameter of order 1. 
There are 0{8~ 1 ) points on the boundary, and the probability that 3 such arms 
exist from one of them is known — at least in the case of a polygon, which is 
enough for our purposes — to behave like S 2 , see e.g. [17] . 

Another possible reason for the non-existence of 3 arms from the source 
of e to the correct portions of the boundary of ft — S (say) is that one of the 
corresponding arms in Q actually lands very close to either A, B or C: Preventing 
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Figure 6. A typical case contributing to P% s (e) - P^~ s s {e). The 
original domain boundary is represented by a solid line, that of the 
shifted domain bu a dashed-and-dotted line; open (resp. closed) 
arms from the source of e are represented as solid (resp. dotted) 
lines, and the additional open path preventing the closed arm from 
connecting to the boundary of fl$ — S is represented as a dashed 
curve. 

it from touching the relevant part of d(Q — S) requires only one additional arm 
from a ^-neighborhood of that vertex — i.e., a total of 2 arms of diameter of order 
1. The probability for that (see [17] also), still in the case when O is a polygon 
with none of A, B or C as a vertex, behaves like 5. Fortunately, there are only 3 
corners on a conformal triangle, so the contribution of these cases is of the same 
order as previously 

Putting everything together, one gets an estimate of the form 



Coming back to our current goal, let £ be the set of oriented edges in a given 
period of fia, and let £ + S be its image by the translation of vector 5. Then, 



with r\ > 0; the existence of such an r\ is ensured by Russo-Seymour- Welsh 
type arguments again, which ensure that, uniformly in e and S, for every edge 



Pts S ( e )= P AAe) [1 + 0(6)]. 



(4.2) 



£ ^(e)PS» = £ 1>{e)P%M [1 + 0(6)} 




= ^^(e)P^( e )+C( ( 5 2+I ') 




\ s (e) = 0(6*>). 



Now, if that is the way the proof starts, what needs to be done is quite clear: 
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• Fix a period E of the graph, 

• Choose a so that the previous sum, over this period, of tp(e)P^ s (e) is equal 
to 0, 

• Use the above estimate to give an upper bound for the same sum on neigh- 
boring periods; 

• Try to somehow propagate the estimate up to the boundary. 

The last part of the plan is the one that does not work directly, because one needs 
of the order of <5 _1 steps to go from £ to 90, and the previous bound is not 
small enough to achieve that; one would need a term of the order of 0(S 3+ri ). It 
is however quite possible that a more careful decomposition of the events would 
lead to additional cancellation, though we were not able to perform it. 

4.2. Using incipient infinite clusters 

Another idea which might have a better chance of working out is based on the idea 
of incipient infinite clusters. We are trying to ensure that ^ip{e)PA,s{e) is equal 
to o{5 2 ) over a period for a suitable choice of a; but for it to be exactly equal to 
depends only on the ratios Pa{s)/ Pa{g') within the period considered, and not on 
their individual values. One can then let 5 go to 0, or equivalently let f2 increase 
to cover the whole space, and look at this ratio. 

Proposition 16. There is a map ir : E(T) — > (0, +oo) such that the following 
happens. Let e, e' be two edges ofT a , which we identify with T for easier notation, 
and let 5 = 1. Then, as increases to cover the whole plane, 



uniformly in the choices of A, B and C on dfl. The map ir is periodic and does 
not depend on the choice of a. 

Proof. The argument is very similar to Kesten's proof of existence of the incipient 
infinite cluster (see [TT]); it is based on Russo-Seymour- Welsh estimates. It will 
appear in an upcoming paper [2]. Notice that there is no requirement for A, B and 
C to remain separated from each other; this is similar to the fact that the incipient 
infinite cluster is also the limit, as n — ► oo, of critical percolation conditioned to 
the event that the origin is connected to the point (n, 0) — which in turn is 
again a consequence of Russo-Seymour- Welsh theory. The speed of convergence is 
certainly different with and without such restrictions on the positions of A, B and 
C, though. □ 

Seeing this Proposition, one is tempted to define a by solving the equation 



where again the sum is taken over one period of the lattice. Indeed, all that remains 
in the sum, over the same period of the lattice, of tp(e)PA(e) is composed of terms 





(4.3) 
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of a smaller order. However, because the limit taken to define n is uniform in the 
choices of A, B and C, in particular it is invariant by re-labelling of the corners 
of the conformal triangle; equivalently, taking Pb instead of Pa leads to the same 
limit. Combining this remark with Equation (|3.2|) . one gets the following identities: 

Ve G E(f), 7r(e) = Tr(r.e) = 7r(r 2 .e). (4.4) 

In other words, n(e) only depends on the source of e. For every edge e = (z, z'), 
let tv(z) := 7r(e): If V is a period of V(T), one has 

by using the remark in Equation l|3.ip . So, the equation l|4.3p is actually always 
true, and does not help in finding the value of a ... 

This is actually good news, because it is the sign of emerging cancellations 
in the scaling limit, which were not at first apparent; that means that the relevant 
terms in (|3.3p are actually smaller than they look at first sight, which in turn 
means that making the leading term equal to by the correct choice of a leads to 
even smaller terms. 

Whether the overall strategy can be made to work actually depends on the 
speed of convergence in the statement of Proposition HH In the case of the tri- 
angular lattice, one can actually use SLE to give an explicit expansion of the 
ratio PA(e) / Pa(s') as fl increases, at least in some cases; this is the subject of an 
upcoming paper [2]. 

5. Other lattices 

5.1. Mixed percolation 

We conclude the speculative part of these notes by some considerations about 
bond-percolation on the planar square lattice. The combinatorial construction we 
perform here does apply to more general cases, but the probabilistic arguments 
which follow do not, so we restrict ourselves to the case of 1? . 

The general idea it to map the problem of bond-percolation on 1? to one of 
site-percolation on a suitable triangulation of the plane. Then, if the arguments 
in the previous section can be made to work, one could potentially prove the 
existence and conformal invariance of a scaling limit of critical percolation on the 
square lattice. 

The key remark was already present in the book of Kesten |10| : For any bond- 
percolation model on a graph, one can construct the so-called covering graph on 
which it corresponds to site-percolation. More specifically, let G\ be a connected 
graph with bounded degree; as usual, let E(G\) be the set of its edges and V{G\) 
be the set of its vertices. We construct a graph Gi as follows: The set V{G%) of 
its vertices is chosen to be E(G\), and we put an edge between two vertices of Gi 
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if, and only if, the corresponding edges of G2 share an endpoint. Notice that even 
if G\ is assumed to be planar, G2 does not have to be — see Figure [7] for the case 
of Z 2 . 




Figure 7. The square lattice Z 2 and its covering graph 



The graph thus obtained from the square lattice is isomorphic to a copy of the 
square lattice where every second face, in a checkerboard disposition, is completed 
into a complete graph with 4 vertices. The next remark is the following: in terms 
of site-percolation, a complete graph with 4 vertices behaves the same way as a 
square with an additional vertex at the center, which is open with probability 
1 — with the same meaning as when we looked at refinement of triangles in 
triangulations, i.e. taking a chain of open vertices in the partially centered square 
lattice and removing from it the vertices which are face centers leads to a chain of 
open vertices in the covering graph of Z 2 . 

So, let again G s be the centered square lattice, as was introduced above, 
and let q £ [0,1]; split the vertices of G s into three classes, to defined a non- 
homogeneous site-percolation model, as follows. Each vertex is either open or 
closed, independently of the others, and: 

• The sites of Z 2 are open with probability p = 1/2; we will call them vertices 
of type I, or p-sites for short, and denote by V\ the set of such vertices; 

• The vertices of coordinates (fc + 1/2, 1 + 1/2) with k + I even are open with 
probability q; we will call them vertices of type II, or g-sites for short, and 
denote by V2 the set of such vertices; 

• The vertices of coordinates {k + 1/2, 1 + 1/2) with k + I odd are open with 
probability 1 — q; we will call them vertices of type III, or (1 — q)-sites for 
short, and denote by V3 the set of such vertices. 

We will refer to that model as mixed percolation with parameters p = 1/2 and q, 
and denote by P\/2, q the associated probability measure. Two cases are of partic- 
ular interest: 

• If q = 1/2, the model is exactly critical site-percolation on the centered square 
lattice G.,; 
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• If g = or g = 1 (the situation is the same in both cases up to a translation), 
from the previous remarks mixed percolation then corresponds to critical 
bond-percolation on the square lattice. 

Besides, all the models obtained for p — 1/2 are critical and satisfy Russo-Seymour- 
Welsh estimates. 

5.2. Model interpolation 

We are now equipped to perform an interpolation between the models at q — and 
q = 1/2. Let (ft, A, B, C, D) be a simply connected subset of Z 2 equipped with 4 
boundary points — say, a rectangle; let U = Uq,,a,b.c.d be the event, under mixed 
percolation with parameters p —1/2 and q, that there is a chain of open vertices of 
fl joining the boundary arcs (AB) and (CD). To estimate the difference between 
the probabilities of U for the two models we are most interested in, simply write 

Pi/2,i/2[U] - Pi/a.oM = J ^ P V2, 3 [^] d 1- (5-1) 

If percolation is indeed universal, then one would expect cancellation to occur, 
hopefully for each value of q; the optimal statement being of the form 

d 

lim sup sup —P 1 / 2 JU] = 0. (5.2) 
otz 2 A,B,c,Dedn <ze(o,i) oq 

The main ingredient in the estimation of the derivative in q is, as one might expect, 
a slight generalization of Russo's formula; to state it, we need a definition: 

Definition 17. Consider mixed percolation on G s , and let E be a cylindrical in- 
creasing event for it (i.e., an event which depends on the state of finitely many 
vertices). Given a realization u> of the model, we say that a vertex v is pivotal for 
the event E if E is realized for the configuration u> v where v is made open, and 
not realized for the configuration lo v where v is made closed. We will denote by 
Piv(E) the (random) set of pivotal vertices for E. 

Proposition 18. With the above notation, one has 

^-Pi/2, q [u] = E 1/2 , q [|Piv(C0 n n n v 2 \ - \pw(u) n n n v 3 \] . 

Proof. The argument is the same as in the proof of the usual formula (in the case 
of homogeneous percolation) ; we refer the reader to the book of Grimmett [8J . □ 



As was the case in the previous section, one can relate the event that a given 
site is pivotal to the presence of disjoint arms in the realization of the model, with 
appropriate color. More precisely, a g-site (say) v £ fl is pivotal if, and only if, the 
following happens: v is at the center of a face of Z 2 ; two opposite vertices of that 
face are connected respectively to the boundary arcs (AB) and (CD) by disjoint 
chains of open vertices; the other two vertices of the face are connected respectively 
to the boundary arcs (BC) and (AD) by disjoint chains of closed vertices; and none 
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of the chains involved contains the vertex v. To state the previous description more 
quickly, there is a 4-arm configuration with alternating colors at vertex v, and the 
endpoints of the arms are appropriately located on dVl — see Figure [H 




Figure 8. A four-arm configuration at vertex v making it pivotal 
for the event U(Q, A, B, C, D). 



The main feature of mixed percolation in the case of the centered square 
lattice is the following: Starting from a configuration sampled according to Pi/2, q 
and shifting the state of all vertices by one lattice mesh to the right, or equivalently 
flipping the state of all vertices, or rotating the whole configuration by an angle of 
7r/2 around a site of type I, one gets a configuration sampled according to P\/2,i- q \ 
on the other hand, rotating the picture by 7r/2 around a vertex of type II or III 
leaves the measure invariant. 

Notice that the existence of 4 arms of alternating colors from a given vertex v 
is invariant by color-swapping; the configuration in Figure[8]is not though, because 
the arms obtained after the color change connect the neighbors of v to the wrong 
parts of the boundary. Nevertheless, one can try to apply the same reasoning as 
in the previous section, as follows: Let v' be the vertex that is one lattice step to 
the right of v. If v is a vertex of type II, then v' is a vertex of type III, and up to 
boundary terms, one can pair all the q-sites in Q to corresponding (1 — g)-sites. 

To estimate the right-hand term in the statement of Proposition [TH1 let 

AO) := P[v G Piv(C/)] - P[v' G Piv([/)]. 

Our goal will be achieved if one is able to show that A(i>) = o( | 1 ~ 1 ) ; or equiv- 
alently, if fls is obtained from a fixed continuous domain by discretization with 
mesh S, if one has 

A(v) = o{5 2 ). 
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In the case of critical site-percolation on the triangular lattice, arguments 
using SLE processes give an estimate to the probability that a vertex is pivotal, 
and from universality conjectures it is natural to expect that they extend to the 
case of mixed percolation on T s . They involve the 4-arm exponent of percolation, 
and would read (still in the case of a fixed domain discretized at mesh S) as 

P[v G Piv(C/)] « (5 5 / 4 . 

So, shifting the domain instead of the point as we did in the last section, one 
would expect an estimate on A(i>) of the order 

A(v) « «5 9 / 4 

(where the addition of 1 in the exponent corresponds to the presence of a 3-arm 
configuration at some point on the boundary on either the original domain or its 
image by the shift). Since 9/4 > 2, that would be enough to conclude. 

However, this approach does not work directly, because of the previous remark 
that the shift by one lattice step does change the measure, replacing q by 1 — q. If 
one is interested in the mere existence of the 4 arms around a vertex, combining 
the shift with color-flipping is enough to cancel the effect; but the estimate one 
obtains that way is of the form 

P[v G Piv{Un,A,B.c.D)} - P[v' G Piv(Uo, B ,c,D,A)] « SP[v G Piv(U Q , A , B ,c,D)} 

(5.3) 

(and Russo-Seymour-Welsh estimates are actually enough to obtain a formal proof 
of this estimate). 

So, once again, what is missing is a way to estimate how much P[v G 
Piv(Us2.A.B,c,D)] depends on the location of A, B, C and D along <9f2; if the 
dependency is very weak, then the estimate in Equation (|5.3p might actually be 
of the right order of magnitude. Once again, it is likely that the way to proceed 
is to use a modified version of the incipient infinite cluster conditioned to have 4 
arms of alternating colors from the boundary, and that the order of magnitude of 
A(v) will be related to the speed of convergence of conditioned percolation to the 
incipient clusters; but we were not able to conclude the proof that way. It would 
seem that this part of the argument is easier to formalize than that of the previous 
section, though, and hopefully a clever reader of these notes will be able to do just 
that . . . 
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